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A. P. Sears, A. Petrenko, G. Catelani, L. Sun, Hanhee Paik, G. 

Kirchmair, L. Frunzio, L. I. Glazman, S. M. Girvin, and R. J. Schoelkopf 

Department of Physics and Applied Physics, Yale University, New Haven, Connecticut 06520, USA 

(Dated: March 7, 2013) 

We study the photon shot noise dephasing of a superconducting transmon qubit in the strong-dispersive limit, 
due to the coupling of the qubit to its readout cavity. As each random arrival or departure of a photon is 
expected to completely dephase the qubit, we can control the rate at which the qubit experiences dephasing 
events by varying in situ the cavity mode population and decay rate. This allows us to verify a pure dephasing 
mechanism that matches theoretical predictions, and in fact explains the increased dephasing seen in recent 
transmon experiments as a function of cryostat temperature. We investigate photon dynamics in this limit and 
observe large increases in coherence times as the cavity is decoupled from the environment. Our experiments 
suggest that the intrinsic coherence of small Josephson junctions, when corrected with a single Hahn echo, is 
greater than several hundred microseconds. 
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Solid-state superconducting quantum systems offer conve- 
nient and powerful platforms for quantum information pro- 
cessing. Rapid progress 0]-[3) is being made in engineer- 
ing qubits and effectively isolating them from the surround- 
ing electromagnetic environment. Despite these efforts, the 
measurement apparatus will always be used to contact the en- 
vironment and is therefore a potential source for decoherence. 

Recently [4] superconducting qubits have been created in- 
side a three-dimensional (3D) resonator, leading to more than 
an order of magnitude increase in coherence time. Inter- 
estingly, the energy relaxation time T\ has increased even 
more than the phase coherence time T 2 *, pointing to a new 
or newly important mechanism for dephasing [5|. These 
devices have a single Josephson junction, eliminating the 
sensitivity to flux noise |6|, and surprisingly show only 
a weak temperature-dependent dephasing, inconsistent with 
some predictions based on extrapolations of junction critical 
current noise (71 ID- In these devices, the qubit state is detected 
by observing the dispersive frequency shift of a resonant cav- 
ity. However, it is known fl9 fflD that in the strong-dispersive 
regime the qubit becomes very sensitive to stray cavity pho- 
tons, which cause dephasing due to their random ac-Stark 
shift lfl2l . It requires increasing care to prevent this extrinsic 
mechanism from becoming the dominant source of dephasing 
as qubit lifetimes increase. Experiments elsewhere |fl3l and 
in our lab fl4| have shown that pure dephasing times can be 
many hundreds of microseconds with careful thermalization 
and more extensive filtering. 

In this Letter, we quantitatively test the dephasing of a qubit 
due to photon shot noise in the strong-dispersive coupling 
limit with a cavity. In this novel regime where the ac-Stark 
shift per photon is many times greater than the qubit linewidth 
7 and the cavity decay rate K (15), the passage of any photon 
through the cavity performs a complete and unintended mea- 
surement of the qubit state. This limit also allows a precise 
determination of the photon number in the cavity using Rabi 
experiments on the photon number-split qubit spectrum (161 . 
With a simulated thermal bath injecting photons into the cav- 



ity and in situ mechanical adjustment of the cavity K, we find 
a pure dephasing of the qubit that quantitatively matches the- 
ory ifTUl . Furthermore, we verify that the qubit is strongly 
coupled to photons in several cavity modes and find that the 
dephasing from these modes accounts for the reduced coher- 
ence times as a function of cryostat temperature. Our mea- 
surements at 10 mK demonstrate that decreasing K leads to 
longer qubit coherence times, suggesting that existing dephas- 
ing in superconducting qubits is due to unintended and pre- 
ventable measurement by excess photons in higher frequency 
modes. 

The experiments were performed (see Fig. [Tji with a trans- 
mon qubit coupled in the strong-dispersive limit to a 3D cav- 
ity, and well approximated by the Hamiltonian (TBI : 

H e ff//j = (O c a*a + ((O q -%aa) b*b- —b f b r bb, (1) 

where the operator a 1 creates a cavity photon and the operator 
b^ creates a qubit excitation. Then (O c is the cavity frequency, 
(Oq and a are the qubit frequency and anharmonicity, and 
X/2n = 7 MHz is the light shift per photon which can be 1000 
times larger than the qubit linewidth of y/27z: = 5 — 12 kHz, 
and the cavity linewidth k/2k = 6 — 120 kHz. The large 
dispersive shift leads to the well-resolved peaks in the qubit 
spectrum shown in Fig. lc, allowing us to conditionally ma- 
nipulate the qubit depending on the cavity photon number 
N |16|. Measuring the height of a given photon number- 
split qubit peak (or the amplitude of a Rabi oscillation at fre- 
quency (O q — Nx) allows a direct determination of the proba- 
bility P(N) for the cavity to have a particular photon number. 

Dephasing of the qubit can be caused by a random change 
in cavity photon number, which shifts the qubit energy by ti% 
per photon and leads to a large rate of phase accumulation 
relative to y. Then the pure dephasing rate y^, obtained in 
a Ramsey experiment for the qubit, depends on the stability 
of the N photon cavity state. When the cavity is connected 
to a thermal bath, the probability P{N) follows a system of 
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FIG. 1: (a) Experimental setup. Noise of varying amplitude at the 
cavity transition frequency is sent into the input port of a 3D res- 
onator. The output port of the resonator has a movable coupler which 
varies the output side coupling quality factor Q c from 1.0 x 10 to 
2.5 x 10 7 . (b) Energy level diagram. The qubit has a transition fre- 
quency that is ac-Stark shifted by — % for each photon in the cav- 
ity, (c) Photon number-splitting of the qubit spectrum. We inject 
noise and create a mean number h of photons in the fundamen- 
tal mode. The peaks correspond to N = 0,1,2 photons from right 
to left, with the cavity Q = 1 x 10 6 , and (*) even without applied 
noise we measure a photon occupation in the TEioi mode of the 
cavity to be h ~ 0.02. (d) Cavity population. Rabi experiments per- 
formed on each photon peak N for increasing noise power with cavity 
Q = 2.5 x 10 5 . The signal amplitude gives the probability of finding 
N photons in the cavity. Two linear scaling factors, fit globally, pro- 
vide conversion from homodyne readout voltage 1 17 1 to probability 
(vertical axis), and from attowatts within the cavity bandwidth to h 
(horizontal axis). Error bars represent lrj fluctuations in the \e) state 
readout voltage. The solid lines are a thermal distribution using the 
fit scaling parameters. 



equations |fl9l for the rate of change into and out of the N pho- 
ton state: dP{N) jdt = Jc(n + 1) (N+ 1 )P(N + 1 ) + KfiNP(N - 
1) — r out P(iV), where the cavity decay rate K = 1/t is the in- 
verse of its decay time T, h is the average number of photons, 
and 



K[(h+\)N + h(N+\)} 



(2) 



combines the spontaneous emission of photons with the stim- 
ulated emission due to thermal photons. Then, in the strong- 
dispersive regime (and neglecting other sources of dephasing) 
the dephasing rate becomes /a = r ou t, and the success of an 
experiment that relies on phase predictability of the qubit re- 
quires a constant photon number in the cavity throughout each 
cycle. 

To verify this prediction for quantitatively, we first cali- 
brate our thermal bath and then obtain K with experiments on 
the photon peaks of the qubit. We can determine the cavity 
decay rate K by exciting the cavity with a short coherent pulse 
while measuring the repopulation of the ground state |g,0) 
(i.e. the amplitude of the zero-photon Rabi oscillations) over 
timescale T. Alternatively, exciting the cavity with a wideband 
noise source that covers the cavity co c transition frequency but 
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FIG. 2: Qubit dephasing due to photon noise, (a) Qubit coherence 
time, determined from Ramsey experiments on the N = or N = 1 
(A) photon peaks, as a function of both cavity Q and n. The dashed 
lines are theory, with an offset due to residual dephasing. Each has 
a slope proportional to K (or 3k for N = 1 experiments), accord- 
ing to Eq. [2] The (o) are coherence times vs. population in TE103 
mode, which also dephases the qubit. (b) Ramsey with no noise in- 
jected, fundamental mode Q = 1 x 10 6 , and r 2 * = 26 /is. The solid 
line is a fit with an exponentially decaying sine, (c) A Ramsey with 
moderate noise. Contrast and r 2 * are reduced. Fundamental mode 



Q = 2.5 x 10 3 



: 0.25, r 2 * = 7.7 [is. (d) Ramsey with high noise. 



Fundamental mode Q = 1 x 10 , n = 3.1, r 2 * = 5.2 fls. Our selec- 
tive (N = 0) pulses produce a loss of contrast and a non-oscillating 
signal addition (orange) as photon population returns to a thermal 
distribution. The dashed black line is a numerical simulation (see the 
Supplemental Materials). 



not the qubit co q transition frequency, creates an average pho- 
ton number h = APbe(T)Q/Q c . Here, A is the linear power 
loss from additional cold attenuation, Pbe = 1 /{e tia l kT — 1) is 
the Bose-Einstein population of the 50 Q. load of the noise 
source at effective temperature T, located outside the cav- 
ity. The total cavity quality factor Q = (OT has an inverse 
which is the sum of the inverses of the coupling quality fac- 
tor Q c of the noise source port, all other port couplings, and 
the internal quality factor Qint- In steady state and for un- 
corrected noise, the probability P(N) of finding the qubit 
in an environment with N photons is a thermal distribution 
P(N) = h N /(n + l)( N+1 \ as verified by the data in Fig. Qi. 
With these measurements we obtain the scaling of n as a func- 
tion of applied noise power for each different value of K\ al- 
lowing a comparison with Eq. 2 using no adjustable parame- 
ters. 

To observe the influence of photon dephasing on our qubit, 
we test Eq. [2] over a wide range of values for both h and K as 
shown in Fig. [2] The photon number is varied by adjusting the 
attenuation following our noise source, while K is controlled 
by retracting the resonator output coupler using a Kevlar string 
connected to the top of the fridge, exponentially increasing 
the Q c as it is withdrawn. For large K, photons enter and 
leave quickly, so long periods uninterrupted by a transit are 
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rare even if the average occupation is low, and the phase co- 
herence time is short. In the Ramsey data shown in Fig.[2]the 
dephasing rate is universally proportional to injected h and K, 
with an offset due to spontaneous decay (if N > 0), and resid- 
ual photons or other intrinsic dephasing. These experiments 
confirm our understanding of the qubit dephasing rate in the 
strong-dispersive limit, and point to the importance of excess 
photons or an effective temperature of a mode for qubit coher- 
ence. 

Importantly, we use slow Gaussian pulses to control the 
qubit in order to exploit the photon-dependence of our Hamil- 
tonian. With a width of (7 = 100 ns, the narrow frequency 
span of the pulses means that Ramsey experiments add sig- 
nal contrast only when the chosen photon number N has re- 
mained in the cavity throughout the experiment, a type of 
post-selection evident in the different scalings of Fig. [2]>d. 
Once conditioned, photon transitions during the experiment 
lead to an incoherent response in our qubit readout, when 
at a random point in time fo an initially prepared superpo- 
sition changes: |y/(f )) = 1/V2(|g,0) + |e,0)) -> |y/(f)) = 
l/\/2(|g,l) +exp[ix(t — to)] |e,l)) for time t > f . Our qubit 
readout [ 17] traces over all photon states, and the unknown fi- 
nal phase of the superposition produces a decay in the Ramsey 
fringes, as the experiment records the qubit excitation despite 
any cavity transition. Additionally, a characteristic bump and 
slope are visible in the data and must be removed before fit- 
ting the Ramsey signal with the usual decaying sine function. 
These features can be understood as the re-equilibration of the 
cavity photon number after the first qubit manipulation condi- 
tionally prepares a certain photon number, and are well fit (see 
Fig. 3 of the Supplement) by a simple master equation which 
includes the incoherent cavity drive as well as qubit and cavity 
decay. 

While the fundamental TEioi mode of our 3D resonator 
serves both as the qubit readout channel and as a mechanism 
for dephasing, the rectangular cavity in fact supports a set of 
TEion modes ll20l whose influence we must consider. Then 
a more comprehensive Hamiltonian than Eq. 1 must incorpo- 
rate many different cavity frequencies, each with a coupling 
strength that depends on antenna length and the positioning 
of the qubit in the cavity [18]. This coupling g„ is large for 
odd-71 TEion modes where the electric field has an antinode 
at the qubit, while the even-n modes have greatly diminished 
coupling to the qubit due to a node along the qubit antenna. 
For our parameters, the fundamental TEioi mode (0\/2n = 
8.01 GHz, (o q /2n = 6.65 GHz, and g x /2n = 127 MHz, the 
qubit anharmonicity a = 340 MHz leads to an ac-Stark shift 
of Xxjln = 7 MHz. Similarly, the first odd harmonics TE103 
with Oh, = 12.8 GHz has a large Xi/ln = 1 MHz. In fact, 
with this mode we can perform high fidelity readout, measure 
the photon mode population (using longer a — 800 ns width 
pulses), and observe its influence on decoherence by inject- 
ing noise near 0)3. In general, we should consider all cavity 
modes that have a non-zero coupling to the qubit as sources of 
significant decoherence. For example, the odd-n TEio„ modes 
at frequency co„ and detuning A„ = co n — a> q , have a coupling 
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FIG. 3: (a) Decoherence due to thermal photons. The coherence 
times extracted from Ramsey (r 2 *) and Hahn echo (^e) experiments 
measured as a function of cryostat temperature. To model dephas- 
ing (dashed lines), we predict population in the TE101 and TE103 
modes of the cavity. Then, we sum the total dephasing rate using 
the measured quality factors for each mode (High Q: T101 = 20 fls, 
T103 = 4 fls; Low Q: %\q\ = 2 [is, Z103 = 400 ns). For high Q, the use 
of a Hahn echo pulse leads to a large Tie because either the photon 
state has much longer correlation time or the remaining dephasing 
similarly occurs at low frequencies. Although the decline in T\ (not 
shown) 1 21 1 contributes to the trend, population in both TE101 and 
TE103 are needed for a good fit. (b) Bose-Einstein population of the 
first two odd-« TE10,, modes at 8 and 12.8 GHz (green) and the coher- 
ence limits they impose individually (blue) and collectively (dashed 
red) for the low Q values measured above. 

gn x \f^n an d an ac-Stark shift %n = 8n a /^n(^n + Ot) which 
decrease only slowly as 1 / n. Consequently, there may be 
many modes with significant dispersive shifts that can act as 
sources of extrinisic qubit decoherence. Moreover, since the 
coupling quality factors of these modes typically decreases 
with frequency, even very small photon occupancies (which 
are usually ignored, not measured or as carefully filtered) must 
be suppressed to obtain maximum coherence. 

The photon shot noise from multiple cavity modes pro- 
vides a simple explanation for the anomalous qubit dephas- 
ing previously observed [4| as a function of cryostat tempera- 
ture. In this case, each cavity mode should be populated with 
the Bose-Einstein probability Pbe and these thermal photons 
can make an unintended measurement of the qubit, disrupt- 
ing phase-sensitive experiments. The predicted occupancies 
for the TE101 and TE103 modes are shown (green lines) in 
the inset of Fig. 3, along with their predicted dephasing (blue 
lines). Having confirmed the dephasing rates for all modes in- 
dividually we can now combine the effect of all modes that 
strongly couple to the qubit: 7^ = £n,K",-. This total ther- 
mal decoherence rate is shown as the red dashed line in the 
inset of Fig. 3, for typical parameters. Since these modes 
have hco„ ^ kgT, the predicted dephasing time is in excess of 
100 microseconds below 80 mK due to the exponentially sup- 
pressed number of blackbody photons. However, since any 
particular mode coupling to the qubit in the strong-dispersive 
limit may have a relatively fast decay time T, even very small 
non-thermal populations n could easily sat- 
isfy hK ^ 1 /27i, limiting the coherence through pure dephas- 
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FIG. 4: Coherence times versus TEini mode decay T. The TE103 
cavity, which naturally decays more strongly through the couplers, 
increases in Q as the entire resonator is decoupled from our coaxial 
lines. While T\ is nearly constant due to the large qubit detuning from 
the cavity, its T£ and Tje increase as the coupling pin is withdrawn 
from the 3D resonator. This is consistent with diminishing dephas- 
ing from cavity modes with K < %, where a photon transit strongly 
measures II II the qubit state. 

ing alone to J 2 * ~ 1 /Yij> =t/n. The measured coherence times 
as a function of temperature are well fit (see Fig. 3) by the 
combined dephasing of thermal occupancy of the TE101 and 
TE103 modes, plus a parameter adjusted to represent the resid- 
ual dephasing in each experiment. This excess could be due to 
another mechanism intrinsic to the qubit, or simply due to in- 
sufficient filtering or thermalization of the apparatus, leading 
to a small non-thermal photon population. 

Further evidence that the true intrinsic coherence limits of 
the 3D transmons have not yet been observed is provided by 
the data shown in Fig. 4, where the qubit relaxation time (T\), 
Ramsey time (r 2 *), and Hahn echo time (T2e) at 10 mK are 
shown as a function of the TE101 cavity decay time. We see 
that the relaxation time is relatively unaffected by cavity life- 
time, since this qubit is sufficiently detuned from the cavity 
to minimize the Purcell effect (Tj. However, we observe a 
general trend where r 2 * and Tie increase as the cavity life- 
time increases. This is consistent with a decoherence due to 
residual photons with ever slower dynamics, but not expected 
due to e.g. junction critical current noise, which should be 
independent of cavity properties. The maximum echo time 
(45 fis w 1.5 Ti) observed here indicates that coherence of 
small Josephson junction qubits is in excess of several hun- 
dred microseconds when corrected by a single Hahn echo. 

In conclusion, we have performed experiments involving 
precise thermal photon populations to quantitatively induce 
qubit dephasing in good agreement with simple theory. We 
find that photons in the fundamental and at least one harmonic 
mode of the cavity strongly couple to a transmon qubit and 
note that at the nominal base temperature of our cryostat they 
produce a negligible amount of dephasing. However, the sen- 
sitivity of the qubit to photons at many frequencies requires 
that we either keep all modes of the cavity in their ground 
state, or else minimize the influence of non-thermal popula- 



tions by reducing their measurement rate ll22l . Inclusion of 
the cavity harmonics in dephasing calculations leads to an un- 
derstanding of the earlier, anomalous, temperature-dependent 
decoherence in our devices |4|. Finally, we find evidence that 
interactions with the residual photons in our 3D cavity likely 
mask the intrinsic coherence time of the Josephson junction, 
whose limits are much longer than qubit coherence times seen 
so far. As qubit linewidths shrink in the future, other effects 
such as quasiparticle parity ||23Tl25l or nuclear spins |26ll may 
further split the qubit spectrum, enabling probes of their state 
dynamics using these procedures. 

We thank Michel Devoret for valuable discussions. L. F. 
acknowledges partial support from CNR-Istituto di Cibernet- 
ica. This research was funded by the Office of the Director 
of National Intelligence (ODNI), Intelligence Advanced Re- 
search Projects Activity (IARPA), through the Army Research 
Office, as well as by the National Science Foundation (NSF 
DMR-1004406). All statements of fact, opinion or conclu- 
sions contained herein are those of the authors and should not 
be construed as representing the official views or policies of 
IARPA, or the U.S. Government. 



[1] A. Houck, J. Schreier, B. Johnson, J. Chow, J. Koch, J. Gam- 

betta, D. Schuster, L. Frunzio, M. Devoret, S. Girvin, et al„ 

Phys. Rev. Lett. 101, 080502 (2008). 
[2] M. Neeley, M. Ansmann, R. Bialczak, M. Hofheinz, N. Katz, 

E. Lucero, A. O'Connell, H. Wang, A. N. Cleland, and J. 

M. Martinis, Phys. Rev. B 77, 180508(R) (2008). 
[3] M. Steffen, F. Brito, D. DiVincenzo, S. Kumar, and M. Ketchen, 

New J. Phys. 11, 033030 (2009). 
[4] H. Paik, D. Schuster, L. Bishop, G. Kirchmair, G. Catelani, 

A. Sears, B. Johnson, M. Reagor, L. Frunzio, L. Glazman, et al„ 

Phys. Rev. Lett. 107, 240501 (2011). 
[5] A. A. Houck, J. Koch, M. H. Devoret, S. M. Girvin, and R. J. 

Schoelkopf, Quantum Inf. Process. 8, 105 (2009). 
[6] F. C. Wellstood, C. Urbina, and J. Clarke, Appl. Phys. Lett. 50, 

772 (1987). 

[7] D. Van Harlingen, T. Robertson, B. Plourde, P. Reichardt, 
T. Crane, and J. Clarke, Phys. Rev. B 70, 064517 (2004). 

[8] J. Eroms, L. C. van Schaarenburg, E. F. C. Driessen, J. H. 
Plantenberg, C. M. Huizinga, R. N. Schouten, A. H. Ver- 
bruggen, C. J. P. M. Harmans, and J. E. Mooij, Appl. Phys. 
Lett. 89, 122516 (2006). 

[9] P. Bertet, I. Chiorescu, G. Burkard, K. Semba, C. Harmans, 
D. DiVincenzo, and J. Mooij, Phys. Rev. Lett. 95, 257002 
(2005). 

[10] J. Gambetta, A. Blais, D. Schuster, A. Wallraff, L. Frunzio, 
J. Majer, M. Devoret, S. Girvin, and R. Schoelkopf, Phys. Rev. 
A 74, 042318 (2006). 

[11] I. Serban, E. Solano, and F. K. Wilhelm, Europhysics Letters 
(EPL)80, 40011 (2007). 

[12] D. Schuster, A. Wallraff, A. Blais, L. Frunzio, R. Huang, J. Ma- 
jer, S. Girvin, and R. Schoelkopf, Phys. Rev. Lett. 94, 123602 
(2005). 

[13] C. Rigetti, S. Poletto, J. Gambetta, B. Plourde, J. Chow, A. Cor- 
coles, J. Smolin, S. Merkel, J. Rozen, G. Keefe, et al., Arxiv 
preprint arXiv: 1202.5533 (2012). 



5 



[14] L. Sun and A. Petrenko, in prep. 

[15] D. I. Schuster, A. A. Houck, J. A. Schreier, A. Wallraff, J. M. 
Gambetta, A. Blais, L. Frunzio, J. Majer, B. Johnson, M. H. 
Devoret, et al., Nature 445, 515 (2007). 

[16] B. R. Johnson, M. D. Reed, A. A. Houck, D. I. Schuster, L. S. 
Bishop, E. Ginossar, J. M. Gambetta, L. DiCarlo, L. Frunzio, 
S. M. Girvin, et al., Nature Phys. 6, 663 (2010). 

[17] M. Reed, L. DiCarlo, B. Johnson, L. Sun, D. Schuster, L. Frun- 
zio, and R. Schoelkopf, Phys. Rev. Lett. 105, 173601 (2010). 

[18] S. Nigg, H. Paik, B. Vlastakis, G. Kirchmair, S. Shankar, 
L. Frunzio, M. Devoret, R. Schoelkopf, and S. Girvin, Arxiv 
preprint: 1204.0587 (2012), Phys. Rev. Lett, in press. 

[19] D. F. Walls and G. J. Milburn, Quantum optics (Springer, 
Berlin; New York, 1994), ISBN 3540588310. 

[20] C. P. Poole, Electron spin resonance : a comprehensive trea- 
tise on experimental techniques (Wiley, New York [u.a.], 1967), 



ISBN 9780470693865. 
[21] G. Catelani, J. Koch, L. Frunzio, R. Schoelkopf, M. Devoret, 

and L. Glazman, Phys. Rev. Lett. 106, 077002 (201 1). 
[22] M. Hatridge and S. Shankar, in prep. 

[23] J. Schreier, A. Houck, J. Koch, D. Schuster, B. Johnson, 
J. Chow, J. Gambetta, J. Majer, L. Frunzio, M. Devoret, et al., 
Phys. Rev. B 77, 180502(R) (2008). 

[24] L. Sun, L. DiCarlo, M. Reed, G. Catelani, L. S. Bishop, D. I. 
Schuster, B. R. Johnson, G. A. Yang, L. Frunzio, L. Glaz- 
man, et al., Arxiv preprint: 11 12.2621 (2012), Phys. Rev. Lett, 
in press. 

[25] O. Naaman and J. Aumentado, Phys. Rev. B 73, 172504 (2006). 
[26] D. Schuster, A. Sears, E. Ginossar, L. DiCarlo, L. Frunzio, 

J. Morton, H. Wu, G. Briggs, B. Buckley, D. Awschalom, et al., 

Phys. Rev. Lett. 105, 1 4050 1 (20 1 0) . 



6 



Supplemental Material for 
Photon Shot Noise Dephasing in the Strong-Dispersive Limit of Circuit QED 

A. P. Sears, A. Petrenko, G. Catelani, L. Sun, H. Paik, G. Kirchmair, 
L. Frunzio, L. I. Glazman, S. M. Girvin, and R. J. Schoelkopf 

Department of Physics and Applied Physics, Yale University, New Haven, Connecticut 06520, USA 

(Dated: June 6, 2012) 



I. EXPERIMENTAL APPARATUS 

The experiment involved a two port 3D resonator with qubit J3, previously reported [1] to have T[ = 42 /is, and 
T,* = T2E = 12 ,us with a fundamental cavity quality factor Q ~ 3.2 x 10 5 . We placed the aluminum 3D resonator at the 
10 mK base plate of a dilution cryostat, attached to a brass translation stage and the copper cold finger of the cryostat. On the 
input side of the resonator, coax lines were attenuated first by 20 dB at the 4 K stage, then at the base plate by a -20 dB directional 
coupler followed by 10 dB of resistive attenuation before the input port to the cavity, similar to our previous experiments (see 
Fig. 1). The output side coupler was attached via a Kevlar string, heatsunk at several stages, to a linear feedthrough at the top of 
the cryostat. With this we were able to translate the coupler stage by up to 3 mm. Hand-formable UT-85 coax cables allow the 
output coupler to be pulled out of the cavity, while a BeCu spring provides a restoring force. This varies the coupling quality 
factor Q c of the port from 1 x 10 5 — 3 x 10 7 , while the input coupler is fixed with coupling quality factor Q m = 1 x 10 7 . After 
the output coupler but still at the base plate, two non-reciprocal isolators attenuated Nyquist-Johnson noise from the wideband 
HEMT amplifier mounted at 4 K while permitting transmission of the readout signal. At room temperature, the signal was mixed 
down to 10 MHz after further amplification. 

II. RABI EXPERIMENTS ON PHOTON PEAKS 

In order to determine the thermal n occupation of a cavity mode, we perform a set of power Rabi experiments, varying the 
power of a single Gaussian pulse of fixed duration applied to one of the N photon qubit peaks. The qubit readout signal after a 
photon selective qubit 7i-pulse is proportional to the difference in qubit excited and ground state population P(e,N) — P(g,N). 
In Fig. 2 each graph contains a power Rabi dataset for N = 0, 1 , and 2, for a different mean number of photons h injected into the 
cavity. The oscillating signal amplitude in a trace is a lock-in measurement of the probability for the cavity to be occupied by N 
photons. Each power Rabi measurement provides a data point in Fig. Id of the main text, with axes for h (horizontal) and P (N) 
(vertical) that have been scaled by the optimal two parameters which fit the whole data set to a thermal distribution. Repeated 
datasets analyzed by the same procedure yielded an h scaling parameter that varied by ±5%. 

Additionally, we found P(e,0) = 5± 1% throughout the experiment, a phenomenon which is not unique [2] among 3D 
superconducting qubits but also not yet completely understood. To obtain this value we compare the contrast of power Rabi 
experiments on the excited qubit peak at m q — a, with and without an initial n inversion on the qubit, similar to measurements 
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FIG. 1 : Block diagram of the measurement setup. 

of cavity photon population. Since our noise source covers even the cavity transitions that have shifted due to the excited qubit 
state (the state dependent shift), this is not expected to interfere with our quantitative measurements of dephasing. We also find 
(see Fig. 2*) even without applied noise we measure a photon occupation in the TEjoi mode of the cavity to be h ~ 0.02. 



III. RAMSEY SIGNAL 



In this section we give a brief explanation of the procedure used to fit the Ramsey experiments. The energy level diagram in 
Fig. lb of the main text can be interpreted as showing that for a given qubit state the cavity is effectively a harmonic oscillator 



8 



0.36 



_ 0.25 
> 

£ 0.20 







0.0* 







0.95 
0.50 ' 



0.0S 











AAAA 



AAAA 



-271 2tc 
Rotation Angle (rad) 



0.1S 



0.25 



0.45 






FIG. 2: Rabi experiments on the TEjni photon number-split qubit. For each curve the qubit Rabi pulse is on a different peak, from N — (red) 
to N — 1 (green) and N — 2 (blue). The vertical scale is the integrated readout signal, and in all graphs P(0) > P{1) > P(2). All graphs have 
the same axes, and the cavity Q is 2.5 x 10 5 . 

with frequency (O c , or co c — % when the qubit is excited. Assuming long qubit relaxation time T\ S> 1/jc, the qubit does not relax 
during the evolution of the cavity occupation. Then for a given qubit state the dynamics of the probability P(N,t) of having N 
photons in the cavity at time r is governed by the rate equations: 



= K(n+l)(N+l)P(N+l) + KnNP(N- 1) - K[n(N + 1) + (n + 1)JV] P(N) 
and the steady-state probability P S (N) is 



1 



(1) 



(2) 



1+n \l+n 

With the qubit initially in the ground state, a jt/2 pulse at the N = peak will coherently split the occupation probability 
between ground and excited states. This means that the initial probability distribution after the jt/2 pulse is 



P(N,0) ■■ 



2 l+n ' 



o. 



N>0 



(3) 
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FIG. 3: A set of Ramsey experiments similar to those used for series in Fig. 2 of the main text. Each graph has the same axes and is labelled 
by the n e jj, as calibrated by the Rabi thermometry experiments. The adjusted data (red) has had the incoherent readout response (orange) 
subtracted from the readout signal, before being fit with a exponentially decay sine (dashed). Additionally, the original data after subtracting 
the final fitted decaying sine is shown (green). The cavity Q is 2.5 x 10 5 . 



After the second Tt/2 pulse at time tf, the probability distribution for N > is simply P(N,tj) as determined by the evolution of 
P(N,t) from the above initial conditions, while P(0,tf) is the usual Ramsey fringe signal [3] which we denote Si{(tf), exhibiting 
decay times which scale according to 1/nK (see Fig. 2 of the main text). As the readout sums over all photon number states, the 
signal S(tf) is 

S(t f )=S R (tf)+P B (t F ) (4) 

with the "bump" Pg defined as 

M'f) = I P(N,t f ). (5) 
N>0 

To fit the signal, we need to calculate the function This can be done by solving the rate equations (1) with the initial 
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FIG. 4: Universal Dephasing Plot: an alternative presentation of our measurements of qubit dephasing versus noise. Here the horizontal axis is 
nK, scaled for each trace, which according to Eq. 2 from the main text we expect to be equal to 7^ for the qubit. The vertical offsets correspond 
to residual dephasing even for n — 0, evident in Fig. 4 in the main text. 



conditions (3). An analytical solution will be given elsewhere. Here we simply give the final answer, 

1 h l — e~ 

2(l+«) 2 l-if^ 

So far we have ignored the qubit and its relaxation. A more detailed analysis gives a reduction of the amplitude of the bump 
(which we account for with the factor M in Eq. 8 below) and its relaxation over the T\ time scale. Then the actual signal is given 
by 

S(t f ) = S R (t f ) + S B (t f ) (7) 

with 

S B {tf)=MP B {t f )e-l T \ (8) 

consistent with the dashed curves in Fig. 3 and in the main text. 

We account for this signal with the following process. First we simulate the system using values obtained via experiment: 
cavity decay time T, and n obtained through calibration. Then, we take the expected Ss(f/) (the excited qubit population with 
N > 0) and remove it from the Ramsey signal, fitting the remainder to a decaying exponential. This typically leads to a correction 
to y that is ~ 10 — 30%. The products of this process are shown for a series with the same K = 1 /t as that in Fig. 2 here, and 
Fig. Id of the main text. 

IV. COUPLED BATHS 



When the resistor component of an isolated RLC circuit is at finite temperature, the harmonic oscillator mode has non-zero 
population h = Pbe(T), the Bose-Einstein occupation at the mode frequency. If other baths are coupled with the oscillator, 



11 



6 



C K XdB n, 

HI — pwvh ^wsHi' 

Att. Z 

4= ^ n eff 



FIG. 5: A harmonic oscillator is coupled via a port to a thermal load, which populates the mode with h. 

then it experiences a total rate of excitation: I* = £n,fQ where ic, for our system includes both input and output pin couplers 
and an additional port representing internal dissipation, while the decay rate of any population in the oscillator is K = £fC;. 
In Fig. 5 we present this as a capacitive coupling C K to a resistor, where we can approximate [4] = m (C K (o) 2 ZoZj nt for 
a resonator with characteristic impedance Z; nt . When a cold attenuator with dissipation X in dB is inserted between a port 
and its thermal load, this effectively reduces Pbe(T) by A = 10 _x / 10 , the linear power reduction amount [5] of the attenuator. 
Then in the steady state £n/K/ = « e ff"' I or « e ff = £re;Ki/£Kj, where n e ff is the mean photon number due to thermal or injected 
noise. For a single external bath at temperature T, the only non-zero n, is injected through the port with decay rate K c , giving 
n eff = APbe(T)Kc/k = APbe(T)QIQ l as in the main text. 
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